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The effects of kinematic cuts on electric charge fluctuations in a gas of charged particles are dis¬ 
cussed. We consider a very transparent example of an ideal pion gas with quantum statistics, which 
can be viewed as a multi-component gas of Boltzmann particles with different charges, masses, and 
degeneracies. Cumulants of net-electric charge fluctuations Xn are calculated in a static and expand¬ 
ing medium with flow parameters adjusted to the experimental data. We show that the transverse 
momentum cut, Pt min < pt < pt max , weakens the effects of Bose statistics, i.e. contributions of 
effectively multi-charged states to higher order moments. Consequently, cuts in pt modify the ex¬ 
perimentally measured cumulants and their ratios. We discuss the influence of kinematic cuts on 
the ratio of mean and variance of electric charge fluctuations in a hadron resonance gas, in the light 
of recent data of the STAR and PHENIX Collaborations. We find that the different momentum 
cuts of Pt min = 0.2 GeV (STAR) and Pt min = 0.3 GeV (PHENIX) are responsible for more than 
30% of the difference between these two data sets. We argue that the pt cuts imposed on charged 
particles will influence the normalized kurtosis kqoq = X 4 /xi °f the electric charge fluctuations. 
In particular, the reduction of k,qoq with increasing pt min will lead to differences between PHENIX 
and STAR data of 0(6%) which currently are buried under large statistical and systematic errors. 
We furthermore introduce the relation between momentum cut-off and finite volume effects, which 
is of relevance for the comparison between experimental data and lattice QCD calculations. 

PACS numbers: 25.75.Nq, 25.75.Gz 


I. INTRODUCTION 


Fluctuations of conserved charges provide a unique op¬ 
portunity to describe the state of matter created in heavy 
ion collisions [J-Q. In particular, higher order cumulants 
of the event-by-event charge multiplicity distribution are 
interesting quantities as they become increasingly sen¬ 
sitive to critical behavior near the phase boundary of 

qcd 3 , 130 . In the Beam Energy Scan (BES) pro¬ 
gram at the Relativistic Heavy Ion Collider (RHIC ) up 
to fourth order cumulants of net-proton fluctuations jl5| , 
as a proxy for net-baryon number fluctuations, and net- 
electric charge fluctuations have been measured 0113 . 

While cumulants of conserved charges exhibit singular 
behavior at the critical point 3 , 0 ] , they are also sensi¬ 
tive probes for pseudo-critical behavior in the vicinity of 
the chiral crossover transition 3 , 10111 , 0 . Further¬ 
more, fluctuations of conserved charges carry information 
on freeze-out conditions in heavy ion collisions [0 0 - l24j 
and they can be influenced by different final state effects 


The proximity of the chemical freeze-out, determined 
by a fit of statistical hadronization models to data 
for particle yields in heavy ion collisions [27], and the 
crossover transition determined in lattice QCD, suggests 
that one may gain access to critical phenomena at the 
QCD phase boundary by measuring fluctuations of con¬ 
served charges at the LHC and in the BES at RHIC 
00 - 

In order to identify critical properties of different fluc¬ 
tuations, one needs to understand a non-critical base¬ 
line for their probability distribution and resulting cu¬ 
mulants. Note that the often used “Poisson baseline”, 
i.e. the Skcllam distribution, is not an adequate ref¬ 
erence for the analysis of net-electric charge and net- 
strangeness fluctuations, due to contributions from multi- 
charged states and quantum statistics effects. This is not 
the case for net-baryon number fluctuations where due 
to the absence of multi-charged baryons and the large 
nucleon mass a Skellam distribution is indeed an appro¬ 
priate reference 0]. 

Net-electric charge fluctuations are dominatecf] by pi- 
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1 This is strictly speaking correct for quadratic and quartic charge 

fluctuations. We will argue later that baryons and kaons also give 

significant contributions to odd order cumulants, e.g. mean and 

skewness. 
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ons. They are therefore strongly influenced by quan¬ 
tum statistics effects. Furthermore, experimental accep¬ 
tance cuts can severely influence electric charge fluctua¬ 
tions. Since the Skellam distribution describes the dis¬ 
tribution for the difference of two independent quantities 
obeying Poisson statistics, limiting the momentum ac¬ 
ceptance does not change the nature of the distribution. 
Thus, fluctuations of net baryon number should be rather 
weakly influenced by kinematic cuts E3- However, due 
to contributions from multi-charged states and quantum 
statistic effects, this is not the case for electric charge 
fluctuations. 

In this study, we discuss the effects of different kine¬ 
matic cuts on the net-electric charge fluctuations in a 
static and expanding medium. We consider the effect of 
a transverse momentum cut-off and of a limited pseudo¬ 
rapidity acceptance on different ratios of the n —th order 
cummulants \n of net-electric charge fluctuations. Sim¬ 
ilar studies, have been carried out in Ref. j29| within a 
hadron resonance gas model. 

In the following, we will focus on the theoretical inter¬ 
pretation of these results, based on a very transparent 
example of an ideal gas of charged pions with quantum 
statistics, which can be viewed as a multi-component gas 
of particles with different charges, masses and degenera¬ 
cies. Such a system contains all relevant properties of 
charge fluctuations that allow for a transparent imple¬ 
mentation and interpretation of different kinematic cuts 
and their consequences. 

In particular, we show that the substantial modifica¬ 
tion of the higher order cumulant ratios with momentum 
cuts , in the acceptance windows of STAR Collaboration, 
that have been discussed in Ref. [2!| may be viewed as 
arising from the suppression of multi-charged states, re¬ 
flecting the Bose momentum distribution in a pion gas. 
We will extend this discussion to different low and high 
transverse momentum cuts, as well as pseudo-rapidity 
windows. We will also focus on properties of the x?/x 2 
ratio which is an important baseline to identify cut-off 
effects on observable which is non-critical with respect to 
the 0(4) universality class. 

We perform a quantitative comparison of our analysis 
with data obtained by the PHENIX and STAR Collabo¬ 
rations. However, in this case we also include the contri¬ 
bution of other hadronic states using a hadron resonance 
gas (HRG) model. We will show, that the increase of 
the lower momentum cut from p tmin = 0.2 GeV (STAR) 

Pimm = 0-3 GeV (PHENIX) does account for a large 
fraction of the differences present in the data for mean 
over variance, Mq/<7q. Similar differences are expected 
to show up in data for the normalized kurtosis, kqctq. 

Finally we will discuss the relation between non-zero 
momentum cut-off and finite volume effects in a non¬ 
interacting pion gas. This allows to argue that lat¬ 
tice QCD calculations of electric charge fluctuations per¬ 
formed in a finite volume are actually adequate to de¬ 
scribe charge fluctuations determined experimentally and 
infinite volume extrapolation may even not be needed 


due to the similarity between finite volume and non-zero 
transverse momentum cut-off effects. 

The paper is organized as follows: In the next section 
we discuss the momentum cut-off dependence of electric 
charge fluctuations in a pion gas. In Section III we intro¬ 
duce STAR and PHENIX data on cumulants of charge 
fluctuations Xn an d their interpretation. In Section IV 
the effects of an expanding medium on fluctuations will 
be discussed. In Section V we compare the effects of fi¬ 
nite size and momentum cuts on Xn ■ Our conclusions 
are presented in Section VI. 


II. CHARGE FLUCTUATIONS IN A PION GAS 

The thermodynamic pressure for a non-interacting 
pion gas at temperature T and electric charge chemical 
potential pq is given by 

p(T,p Q )=-Tj^j ^hr(l-e-^-^/ T ) (1) 


with 


7?max P*max 

J d 3 p = 2n J dr] J dp t p t \p\ (2) 

Vtnin Pt min 

where E p = i/|p| 2 + m 2 is the pion energy and Qi = 
0, ±1 denotes their electric charge. 

The momentum integration in Eq. ©> is expressed 
in terms of the transverse momentum p t and pseudo¬ 
rapidity rj = [ln(|p| +p z )/(\p\ ~ Pz)]/ 2- We consider the 
experimental acceptance window, 0.2 GeV < p t < 2 GeV 
and \r]\ < 0.5 for STAR [lgj, and 0.3 GeV < vt < 2 GeV 
and \ri\ < 0.35 for PHENIX measurements ~fl7| . 

The cumulants of electric charge fluctuations are ob¬ 
tained from 


,q_ 3" p(T, p Q ) 

n d(p Q /T) n T 4 


( 3 ) 


To characterize the properties of the system, it is con¬ 
venient to introduce different ratios of Xn/Xm■ Here¬ 
after we focus on two ratios, the mean over variance, 
Mq/&q = X? 1 X 2 and the normalized kurtosis kqVq = 

Xa/X2- 

In the presence of a non-zero momentum cut-off the 
cumulants may be represented by the standard cluster 
expansion 

2 00 

Xn = k n ~ 2 k 2 {km/T) cosh (kp Q /T), n = even 

k=l 
2 00 

Xn = ^2 k n ~ 2 k 2 {km/T) sinh (kp Q /T), n = odd 
k =1 


( 4 ) 
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FIG. 1. The lower transverse momentum cut-off dependence 
of the first order cumulant x? °f the electric charge fluc¬ 
tuations in a Bose gas at different values of the particle 
mass, m = m n ,mK,m p and rriK * at T = 160 MeV and 
fiQ = — 1 MeV. The results are normalized to the correspond¬ 
ing values at p tmin = 0. 

where the function K 2 (x) is defined as, 

V max Pfmax 

k 2 {km/T) = J dr] J dp t p t \p\e~ k X , M 2+rn2 / T . 

Vtnin Pt min 

( 5 ) 

It reduces to the modified Bessel function K 2 (km/T ) 
after integration over the entire momentum space, —oo < 
r\ < oo and 0 < pt < oo. This makes the influence of 
the Bose statistics very transparent and highlights the 
similarity between cut-off effects on Bose particles and 
multi-charged hadrons. 

From Eq. 0 it is clear that in a Bose gas \n can be 
regarded as a sum over contributions from multi-charged 
particles with mass km, charge k and degeneracy fc ra-4 , 
where k = 1, 2, • • • , oo. The k = 1 term corresponds to 
the Boltzmann approximation for which cumulant ratios 
Xn+ 2 /Xn = 1; be., the normalized kurtosis kq<Jq = 1. For 
k = 1 the corresponding probability distribution P(N) 
for net charge N is the Skellam function mm. Contri¬ 
butions of multi-charged Boltzmann particles with k > 1, 
result in a wider P(N), which has a broader tail than 
the Skellam function |28|. This leads to larger quar- 
tic (xi) than quadratic (x 2 ) charge fluctuations. More 
generally, the higher order terms in k, which become in¬ 
creasingly important for large n owing to the degeneracy 
factor fc™ -4 , imply that Xn+2 > Xn ■ Consequently, in 
the presence of multi-charged particles, as well as parti¬ 
cles with Bose statistics, the normalized kurtosis kqOq 
is always larger than unity in a gas of non-interacting 
hadrons. In Eq. 0 the momentum cut-off dependence 



FIG. 2. The lower transverse momentum cut-off dependence 
of the quadratic [x^, (a)] and quartic [\J, (b)] cumulants of 
the electric charge fluctuations in a pion gas at different values 
of the temperature, 100 MeV < T < 200 MeV and for fiQ = 
— 1 MeV. The results are normalized to the corresponding 
values obtained at Pt min = 0. 

enters through the function K 2 {km/T ) where the effect 
of pt min 0 quantitatively depends on m , T and k sep¬ 
arately. A non-zero pt min thus modifies the relative con¬ 
tribution of different fc-terms to \n- 

Fig- CD displays the pt min dependence of the first cumu¬ 
lants at T = 160 MeV and pq = — 1 MeV for various 
particle masses, to = m v , mx, m p and normalized 

to their values calculated with pt min = 0. 

The charge chemical potential is assumed to be nega¬ 
tive and small, as expected in heavy ion collisions at high 
energies. For m = mx* , the mass is already so large 
that the quantum statistics effect is almost negligible, 
thus the result can be regarded as that of the Boltzmann 
approximation where the higher fc-terms in Eq. 0 are 
neglected. 

By decreasing the particle mass, one sees a stronger 
dependence of xf on Pt min , reflecting the contribution 
from k > 1 terms which are also understood as the en¬ 
hancement of the particle number at low pt in the Bose 
distribution. 
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In Fig. [2] we show the cut-off dependence of the two 
even order cumulants X 2 an d X? for a pion gas at 
Hq = — 1 MeV and at different temperatures. It is clear 
that cut-off effects are more severe for higher order mo¬ 
ments and they are more pronounced at lower tempera¬ 
ture. This temperature dependence can be deduced from 
the K 2 {km/T) function by changing the integration vari¬ 
able from pt to x = kpt/T , which shifts the lower limit of 
the integration range to £min = kp ttairi /T. An increase of 
the temperature, therefore effectively lowers the p t cut. 
Consequently, the cumulants at higher temperature take 
a value closer to those without p t cut. 

The difference in the cut-off dependence of X 2 an d Xa 
will also influence their ratio, Xa IX 2 = K Q a Q- However, 
the effect is weaker in kqctq, due to a partial cancela¬ 
tion of a similar dependence of quartic and quadratic 
fluctuations on p t momentum cut-off. This is evident 
in Fig. [3] where we show the momentum cut-off depen¬ 
dence of X 4 1 X 2 ■ The temperature dependence of kq<Tq 
indicates, that there are stronger cut-off effects at higher 
temperature, where the multi-charged states in the clus¬ 
ter sum become increasingly important. 

The above discussion of the temperature dependence 
of the cumulants and their ratios at different pt cuts can 
be extended to more general cases. For \hq\ <C T and 
\hq\ -C to, which is the case in heavy ion collisions at 
RHIC and LHC energies, the cluster expansion Eq. 0 
implies, that 

X2n+l/X2n ~ Hq/ T > ( 6 ) 

irrespective of the p t cut. This is because, the I< 2 {km/T) 
functions cancel in these ratios. Consequently, in the 
pion gas, the lowest ratio X\ /X 2 does not depend on the 
Pt-cut. Similarly, one finds 

Q Q 

X2n+1 n . X2n+2 rr-?\ 

—Q “ Q ’ 1C 

X2n—1 X 2 n 

which implies that X 4 /X 2 — Xs/Xi- These relations 
hold exactly for hq = 0, but deviations at hq ^ 0, in the 
range of RHIC and LHC energies, are found to be less 
than 1%. 

Fig.|4]shows the normalized kurtosis, kqctq, calculated 
in a pion gas at T = 160 MeV as a function of the pseudo¬ 
rapidity cut ?7 m ax, and for different pt cuts corresponding 
to the STAR and PHENIX experiments, respectively. 

From Figs. EH we conclude that the dependence of 
electric charge cumulants on transverse momentum cuts 
is strong while the dependence on different pseudo¬ 
rapidity cuts is weak. We also checked the dependence 
of cumulant ratios on the high momentum cut-off, , 
and found that this is negligible for pt max > 1 GeV. (see 
Sec.HVl) 

The difference in the dependence of kq<Jq on the trans¬ 
verse momentum and pseudo-rapidity cuts is transparent 
from the decomposition of the Bose integral in terms of 
the multi-charged cluster sum. A transverse momentum 



FIG. 3. The lower transverse momentum cut-off dependence 
of the ratio of quartic and quadratic electric charge fluctu¬ 
ations in a pion gas at different values of the temperature, 
100 MeV < T < 200 MeV for p Q = - 1 MeV. 
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FIG. 4. Normalized kurtosis kq °q = X 4 /X 2 for an id eal 
pion gas at T — 160 MeV, HQ = 0 and in different trans¬ 
verse momentum windows. Open circles stand for the fully 
integrated pt range. Open triangles and squares show the 
result with the same pt acceptance as STAR and PHENIX 
measurements, respectively f~if|IT?]. 


cut-off, _Pt min , increases the effective mass of particles, 
m e ff = ^ jm 2 + Pt min , and thus further suppresses the 
statistical weight of the fc-cluster contribution relative 
to the leading Boltzmann term. Consequently, increas¬ 
ing pt min reduces kqctq towards unity, which is its value 
for the Boltzmann statistics. Cuts in pseudo-rapidity, on 
the other hand, do not restrict the minimal momentum 
in the particle dispersion relation and therefore do not 
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lead to an additional suppression of multi-charged clus¬ 
ters contributing to KqcTq- 


III. ELECTRIC CHARGE FLUCTUATIONS IN 
THE PHENIX AND STAR DATA 

We have shown in the previous section, that a cut-off 
imposed on the low momentum part of the particle phase- 
space can have significant effects on different cumulants 
of net electric charge fluctuations and their ratios. Thus, 
the difference of 100 MeV in the transverse momentum 
cut-off used by the PHENIX and STAR Collaborations 
El 113 , should naturally imply differences in their results 
on cumulants of electric charge fluctuations. 

However, in order to quantify the influence of momen¬ 
tum cut-off effects on charge fluctuations in the context 
of heavy ion data, one needs to extend the pion gas calcu¬ 
lation by including contributions from all known charged 
hadrons and resonances. In the following, we will apply 
a hadron resonance gas (HRG) model which is very suc¬ 
cessful in the description of the equation of state and dif¬ 
ferent fluctuation observables calculated in lattice QCD 
in the hadronic phase [Hj-|34|]. 

We first consider the influence of a non-zero pt min on 
the ratio of mean and variance of net-electric charge fluc¬ 
tuations, Mq/ctq = Xi/X 2 - This quantity is statisti¬ 
cally well under control and a direct comparison between 
data obtained by STAR and PHENIX, as well as with 
model calculations, is meaningful. 

The ratio Mq/oq involves the first order cumulant 

Xi , which receives significant contributions from strange 
mesons as well as charged baryons. In fact, in the case 
relevant for heavy ion collisions, i.e. in a strangeness 
neutral system, Ms = 0, with a fixed ratio of net-electric 
charge and net-baryon number, Mq/Mb — 0.4, the elec¬ 
tric charge chemical potential required to fulfill these con¬ 
straints is negative. The contribution of charged, non- 
strange mesons to Mq thus is negative and the mean 
becomes positive only due to contributions from charged 
baryons and strange mesons. Although in a single com¬ 
ponent Bose gas the ratio Mq/oq does not depend on the 
pt cut-off, Eq. ©, the cumulants themselves are modi¬ 
fied depending on the mass of this component, as seen 
in Fig. [I] However, in a multi-component system, the 
change of Mq and Oq of each hadronic contribution does 
not cancel anymore in their ratio. 

We have calculated Mq/ctq in a HRG model using 
three different scenarios for introducing the low momen¬ 
tum cut-off: (A) pt min is introduced only for the pion 
contribution, (B) an identical cut-off is introduce for 
all charged hadrons, and (C) a larger cut-off, Pi min = 
0.4 GeV, is introduced for protons and anti-protons, as it 
is done in the electric charge measurements of the STAR 
Collaboration. 

We show results from this analysis in Fig. 0 where 
Mq/ctq is calculated along the freeze-out line in central 
heavy ion collisions, with the energy dependent thermal 
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FIG. 5. The ratio of mean and variance, Mq/g'q, of elec¬ 
tric charge fluctuations for the three different scenarios (see 
text) with STAR acceptance (open symbols) and PHENIX 
acceptance (closed symbols) in a HRG model. The energy 
dependence of thermal parameters at chemical freeze-out are 
taken from Ref. E3- The results are normalized with the 
values for p tmin = 0. 


parameters, hb, HQ, Hs and T from Ref. E3- We find 
that the effect of a low momentum cut-off on Mq/oq 
is almost independent of the collision energy. Compar¬ 
ing the different low momentum cut-off scenarios, one 
finds that the set (A) gives the largest effect, since it 
corresponds to cutting off the negative contribution of 
pions to Mq and the pion contribution to erg is reduced. 
Incorporating the cut-off for other hadrons reduces the 
changes from the values without the cut-off. This is be¬ 
cause the kaon contribution partly compensates the in¬ 
crease of Mq due to the change of the pion contribution. 
The result of set (C) shows a weak ^/sNN dependence 
below y/SNN < 20 GeV, due to the dominance of the 
proton contribution. The difference among different sce¬ 
narios is smaller than that observed between STAR and 
PHENIX acceptance. 

We conclude that a non-zero transverse momentum 
cut-off increases the ratio Mq/cfq over the value obtained 
without any cuts. The larger p tmin used by PHENIX 
suggests that their Mq/ctq ratio is about 10% larger 
than the STAR values. This accounts for about 1/3 
of the difference between the PHENIX and STAR data. 
The ratio of these two data sets is, within large er¬ 
rors, constant over the entire energy range. A fit yields 
(Mg/CTg)pHENix/(AfQ/crg)sTAR = 1.35(5). Note, how¬ 
ever, that more than 99% of the error quoted by PHENIX 
is systematic. 

In Fig.[6]we show the data measured by STAR E3 and 
PHENIX []J. They have been normalized by the corre¬ 
sponding ratio for net-proton fluctuations measured by 
STAR [l]|. As can be seen from the right hand panel of 
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FIG. 6. The ratio of mean and variance, Mq/ctq, of electric 
charge fluctuations obtained by STAR (circle) and PHENIX 
(squares) Collaboration at different beam energies Eta, 
(left panel). Data have been normalized to the corresponding 
ratio for net-proton number fluctuations, Mp/ap, obtained 
by STAR 0. The right hand panel shows the rescaled data 
(Mq/oq )p tmin =o = r(AfQ/4)n mia - Where the rescaling fac¬ 
tor r is taken from the result of set (C) for the STAR data 
and set (B) for the PHENIX data shown in Fig. [5] 


Fig. [6j rescaling the STAR and PHENIX data to vanish¬ 
ing transverse momentum cut-off values removes a large 
part of the differences between these data. In fact, within 
the large systematic errors of the PHENIX data it yields 
consistent results. 

In Fig. [7] we show the influence of a non-zero transverse 
momentum cut-off on the normalized kurtosis, kq<Tq. 
Compared to the pion gas results shown in Fig. [3j the 
reduction of kqctq is smaller and turns out to be 10-20% 
due to a smaller cut-off dependence of the heavier addi¬ 
tional degrees of freedom in a HRG. The small discrep¬ 
ancy in the different scenarios nevertheless implies that 
the main reason for the reduction is the pion effect. Con¬ 
trary to Mq/ctq, this result shows a stronger dependence 
of kqOq on y/SNN- Thus, we expect that the difference in 
-Ptmin between STAR and PHENIX acceptance will lead 
to about 6% difference in kq<Jq at ^snn = 200 GeV, 
and it will be negligible at ^ snn = 7.7 GeV. We note, 
that similar results have been obtained in Ref. [29| where 
a further reduction of kq&q due to effects of resonance 
decays is also discussed. 



FIG. 7. Normalized kurtosis kqctq for the three difference 
cut-off scenarios in a HRG model along the freeze-out line. 
The legends are the same as in Fig. [o] 
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FIG. 8. The pt distribution of charged pions. The spectrum of 
the STAR data [35j is shown as triangles. The dashed-line is 
the model results at chemical freeze-out (T = 160 MeV) from 
a longitudinally expanding fireball. The solid-line represents 
the distribution at thermal freeze-out for a longitudinally and 
transversally expanding fireball with parameters fitted to the 
STAR data fH. 


IV. EFFECTS OF EXPANSION 

Until now, we have discussed the influence of momen¬ 
tum cuts imposed in a thermal medium at rest. In 
heavy ion experiments, however, the kinematic cuts on 
measured fluctuations are applied to hadrons that move 
freely only after kinetic freeze-out. Consequently, for the 
comparison of experimental data with theoretical calcu¬ 
lations, it is important to verify how momentum cuts im¬ 
plemented in different stages of the fireball evolution can 


influence the measurement of fluctuations of conserved 
charges [29ij . 

Here we employ an expanding source model with pa¬ 
rameters chosen such that the pion p t spectrum measured 
by the STAR Collaboration [35| is reproduced. We com¬ 
pare the effect of a pt cut imposed either at the time of 
chemical or thermal freeze-out on the measured values 
for the kurtosis of net-electric charge fluctuations. 

Assuming boost invariance along the collision axis and 
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transverse expansion with a linear rapidity profile, as well 
as an instantaneous emission at t = Tf hyper-surface, we 
parametrize the phase space distribution of pions as |.36l | 


S(x,p) 


m T cosh(y ~ r) s )5(T - Tf) -* 2 +f 

-- n B {u-p)e -55*-. 


( 8 ) 

Here y = ln[(£J p + p z )/(E p —p z )}/2 is the rapidity of the 
pions, p s = \n[(t + z)/(t — z)\/2 is the space-time rapidity 
which is equal to the longitudinal flow rapidity yL , and R 
is the Gaussian transverse size of an expanding fireball. 
The flow 4-velocity is parametrized as 


u^ = (cosh yx cosh p s , sinh yx cos </>, sinh yx sin (j>, (9) 

cosh yx sinhjjg). 


The transverse flow rapidity is assumed to be yr = 
ppr/R , with r = \Jx 1 + y 2 and a flow control parameter 
rjf which is adjusted so as to reproduce the p t spectrum 
of pions. 

Since the invariant spectrum is obtained from the 
source function from Eq. m , as 

e p^= J d 4 xS(x,p), ( 10 ) 

the p t spectrum can be obtained by integrating over 
space-time variables, as f drx = 2iiTf J dr J dry, f rdr. 
Consequently, cumulants of charge fluctuations are given 

by 


(il) 

OP'Q J &P 

The proportionality factor is canceled out when taking 
ratios such as Xa/xS- Note that Eq. (fill) reduces to 
Eq. (O when considering a static source, i.e. by choosing 
yL = 0 and yx = 0. 

Fig. [8] displays the pt spectrum of n~ measured by 
the STAR Collaboration (Hsj . We fit the data with the 
above source model with T = 120 MeV and R = 6 fm, 
and obtain pj = 0.373. Although the overall scale factor 
is also fitted, our result still underestimates the yields 
at the lowest pt bin. For comparison we also plot the 
spectrum at T = 160 MeV without transverse flow, i.e. 
for yx — 0 but at finite yL- One sees that it does not 
reproduce the slope of the data. Therefore, it is expected 
that the influence of a p t cut on fluctuation observables 
can be different at chemical and thermal freeze-out due 
to differences in the pt spectrum. 

Fig. [9] shows the normalized kurtosis kqOq calculated 
in an expanding and static source with and without 
transverse flow for different Pt min . The pseudo-rapidity 
window is fixed to |7y| < 0.35 and the upper momentum 
cut to p* = 2 GeV. For p t . < 0.2 GeV there is a 

x <-max r i-min 

clear reduction of Xa/ x? relative to its value obtained 
with a static source from Eq. ©, when the longitudinal 
and/or transverse expansion is included. However, for 
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FIG. 9. The kurtosis kqoq = \a /x? as a function of pt min - 
The long-dashed line stands for the static source from Eq. [T] 
The short-dashed and solid lines denote the results obtained 
from the spectra shown in Fig. [S] 


Pt m i„ St 0-2 GeV the effect of a p t - cut on kqctq is simi¬ 
lar in all cases. The values of the kurtosis do not differ 
much for pt min > 0.2 GeV and the pt min dependence is 
only slightly modified in the presence of transverse flow. 
As can be seen in Fig. [8] the p t spectrum at thermal 
freeze-out, T = 120 MeV and pf = 0.373, clearly has a 
larger slope than that calculated at T = 160 MeV with¬ 
out transverse flow. This difference in the spectral shapes 
leads to a somewhat slower convergence of acqOq to unity 
with increasing pt min • Nonetheless, the resulting decrease 
of kqctq in the interval pt min £ [0.2 GeV, 0.3 GeV] is as 
large as that found for the static source. 

For completeness we also show in Fig. [TO] the influence 
of Pimax on k Q°q at fixed pt mln = 0.3 GeV and \p\ < 0.35. 
There is no notable change of the kurtosis when decreas¬ 
ing Ptmax from 2 GeV to 1 GeV. For p tmax < 1 GeV, the 
relative contribution of lower momentum particles with 
charges k > 1 in Eq. ©, i.e. higher order terms appear¬ 
ing due to Bose statistics, start to dominate fluctuations, 
consequently kqctq is increasing with decreasing pt max . 
The dependence of kq<Jq on pt max is similar for a static 
and transversally expanding medium. 

The above results on flow effects were derived in a 
rather simplified system of a non-interacting pion gas. 
However, due to leading contributions of pions to elec¬ 
tric charge fluctuations and a large mass of multi-charged 
baryons, our conclusions should also hold in a hadron res¬ 
onance gas which accounts for contributions of all known 
hadrons and resonances. 

The sensitivity to kinematic cuts may still be modi¬ 
fied if charge fluctuations are influenced by the critical 
chiral dynamics [37 ], which is clearly not included in this 
study. For the case of charge fluctuations in the crossover 
region, however, we conclude that flow effects can safely 
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FIG. 10. Dependence of the ratio of quartic (x?) and 
quadratic (x?) cumulants of electric charge fluctuations on 
Pt max f° r fixed pt min = 0.3 GeV and |? 7 | = 0.35. The results 
correspond to pion spectra shown in Fig. [8] 

be ignored in the analysis of charge fluctuations. 


V. TRANSVERSE MOMENTUM CUT-OFF 
AND FINITE-SIZE EFFECTS 

In the previous sections we have seen that a quanti¬ 
tative comparison between experimental data taken in 
specific acceptance windows, in particular with a non¬ 
zero transverse momentum cut-off, with the thermal res¬ 
onance gas model requires also to take into account cor¬ 
responding cut-offs in the model calculation. 

In a non-interacting hadron gas a non-zero pt cut re¬ 
duces the available phase space, which may be viewed as 
an increase of the effective mass of particles, and thus 
leads to corresponding changes in thermodynamic ob¬ 
servables. This effect is similar to the thermodynamics 
of particles in a finite volume fdii [il9| . The finite size of a 
system, e.g. a box of size L 3 , suppresses low momentum 
states. 

When considering a box with periodic boundary condi¬ 
tions the lowest non-zero momentum of bosons in such a 
box is 27 x/L. In the case of a non-interacting scalar field 
theory, which describes the thermodynamics of spin zero 
bosons, i.e. pions, it is straightforward to analyze the 
volume dependence of thermodynamic observables [38j |. 
In particular, one can calculate the volume dependence 
of conserved charge fluctuations and compare the influ¬ 
ence of finite volume effects with those of a transverse 
momentum cut-off. 

In Fig.[lT]we show a calculation of quadratic and quar¬ 
tic charge fluctuations of a pion gas in a finite volume. 
Indeed, there is a clear correspondence between finite 
volume and momentum cut-off effects in such a non- 




FIG. 11. Quadratic [xj, (a)] and quartic [x?, (b)] cumu¬ 
lants of electric charge fluctuations in a non-interacting pion 
gas calculated in a finite box of size L 3 for different m n /T 
(symbols). Results are compared to calculations in an infi¬ 
nite box but with a non-zero transverse momentum cut-off 
Pt min = 2nk/L with k = 0.707 (lines). 

interacting gas. As can be seen from Fig.[lTJ all quadratic 
and quartic charge fluctuations in a volume of size L 3 can 
be reproduced by a calculation of charge fluctuations in 
an infinite box, but with a transverse momentum cut-off 

Ptmin = 27 xk/L with k = 0.707 . (12) 

The proportionality factor k is found to be rather weakly 
depend on the pseudorapidity acceptance. Indeed, for 
| A 77 I = 0.35 one gets k = 0.707, while for |A 77 I = 5 one 
finds k ~ 0.5. 

In a non-interacting pion gas the scale for all dimen- 
sionfull quantities is set by the temperature. Conse¬ 
quently, the size of a finite volume or a non-zero mo¬ 
mentum becomes meaningful in units of T. We thus 
may interpret the above results in terms of LT , which 
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is conventionally used in lattice field theory. In par¬ 
ticular, lattice QCD calculations are performed on lat¬ 
tices of size L 1 2 3 4 T _1 = (N a a) 3 N T a and the continuum 
limit ( N a — > oo,N t — > oo) is taken for fixed aspect ra¬ 
tio N a /N T = LT. Typically LT ~ 4 in such calcula¬ 
tions. We may compare this to the momentum cut-off 
p tmin ~ (0.2 — 0.3) GeV currently used in heavy ion ex¬ 
periments. At a freeze-out temperature T ~ 160 MeV 
these momentum cut-offs expressed in units of temper¬ 
ature are, Pt min /T — (1.25 — 1.875). Using Eq. (fl2ll . 
which relates Pt min /T to a box size LT, we find that 
these cut-offs correspond to aspect ratios, LT ~ 2.4 —3.6. 
Thus, lattice calculations performed with an aspect ratio 
LT = 4 seem to be appropriate to characterize the ther¬ 
modynamics probed with transverse momentum cut-offs 
similar to those used by the STAR Collaboration. 

VI. CONCLUDING REMARKS 

We have discussed the influence of kinematic cuts in 
transverse momentum and pseudo-rapidity on electric 
charge fluctuations in a static and expanding medium. 

Our basic arguments were given by considering a very 
transparent example of an ideal gas of charged pions with 
quantum statistics. Such a system can be viewed as a 
multi-component gas of Boltzmann particles with differ¬ 
ent charges, masses, and degeneracy factors, thus it con¬ 
tains all relevant properties to study charge fluctuations, 
with a transparent implementation of different kinematic 
cuts. 

We demonstrate, that a non-vanishing lower p t cut in 
the net-electric charge fluctuations, can substantially re¬ 
duce the second and fourth order cumulants of net elec¬ 
tric charge fluctuations, as well as, their ratio quanti¬ 
fied by the normalized kurtosis kq<Jq. These conclusions 
are valid in a static and expanding medium. The ef¬ 
fect of pseudo-rapidity cuts was shown to be rather neg¬ 
ligible. Consequently, the observed differences between 
STAR and PHENIX data on cumulants of electric charge 
fluctuations could be traced back to their different p t ac¬ 
ceptance. 

The experimental data on the ratio of mean and vari¬ 
ance, Mq/oq, and normalized kurtosis kq<Tq, were ana¬ 
lyzed in a hadron resonance gas. It was shown, that in 
heavy ion collisions, the influence of a non-zero transverse 


momentum cut-off, Pt min , on Mq/uq is almost indepen¬ 
dent of the collision energy, and that it increases, with 
increasing p tmin . Thus, the larger p tmin used by PHENIX, 
suggests that their Mq / ctq should be larger than the cor¬ 
responding STAR values. We have shown, that indeed, 
by rescaling the STAR and PHENIX data to vanish¬ 
ing transverse momentum, makes these data compatible 
within the still large systematic errors of the PHENIX 
data. Furthermore, it was argued that the reduction of 
kqctq due to the shift of pt min will lead to differences 
between PHENIX and STAR data of 0(6%), which cur¬ 
rently are buried under large statistical and systematic 
errors. 

Having in mind a possible comparison between experi¬ 
mental data on fluctuations of conserved charges and lat¬ 
tice QCD results, we have discussed the relation between 
momentum cut-off and finite volume effects. We have 
shown, in the context of an ideal pion gas, that cumu¬ 
lants of electric charge fluctuations in a finite box of size 
L 3 , can be reproduced by a calculation of charge fluctua¬ 
tions in an infinite box, but with a transverse momentum 
cut-off pt min = 27 rk/L, where A: is a constant. 

Finally, considering typical box sizes used in lattice 
QCD calculations, and the lower transverse momentum 
cut-off scale in the experimental data on electric charge 
fluctuations, we have argued, that their direct compari¬ 
son at a chiral crossover temperature, might be justified. 
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